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Abstract 

Phase diagram of tlie spin-1 quantum Heisenberg model with both exchange as weU as single- ion 
anisotropy is constructed within the framework of pair approximation formulated as a variational 
procedure based on the Gibbs-Bogoliubov inequality. In this form adapted variational approach is 
used to obtain the results equivalent with the Oguchi's pair approximation. It is shown that the 
single-ion anisotropy induces a tricritical behaviour in the considered model system and a location 
of tricritical points is found in dependence on the exchange anisotropy strength. 
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I. INTRODUCTION 



Phase transitions and critical phenomena of quantum spin systems currently attract a 
great deal of interest [l]. As usual, the quantum Heisenberg model is used as a basic gener- 
ating model which should be appropriate for investigating quantum properties of insulating 
magnetic materials ^, y, Q, [s]. However, a rigorous proof known as the Mermin- Wagner 
theorem fG] prohibits a spontaneous long-range order for the isotropic spin-1/2 Heisenberg 
model on the one- and two-dimensional lattices and hence, the spontaneous orderingmight 



in principle appear either if a three-dimensional magnetic structure is considered 



3] or a 



non-zero magnetic anisotropy is involved in the studied model Hamiltonian [s^. On the 
other hand, it is currently well established that obvious quantum manifestations usually 
arise from a mutual combination of several factors, especially, when the low-dimensional 
magnetic structure is combined with as low coordination number as possible and low quan- 
tum spin number. Apparently, these opposite trends make hard to find a long-range ordered 
system that simultaneously exhibits evident quantum effects. Investigation of quantum spin 
systems, which can exhibit a non-trivial criticality, thus remains among the most challenging 
tasks in the statistical and solid-state physics. 

Over the last few decades, there has been increasing interest in the study of the effect of 
different anisotropies (single-ion, Dzyaloshinskii-Moriya, exchange) on the critical behaviour 
of the spin-1 quantum Heisenberg ferromagnet. The main interest to study this model 



system arises since Stanley and Kaplan j^, llO| proved the existence of a phase transition 
in the two- and three-dimensional Heisenberg ferromagnets. In addition, the ferromagnetic 
quantum Heisenberg model with the spin-1 has relevant connection with several nickel-based 
coordination compounds, which provide excellent experimental realization of this model 
system jlll. [l^. Up to now, the spin-1 quan tum Heisenberg model has been explored within 
the standard mean- field approximation [isl. 11^. random phase approximation 15 [or linked- 

2l|, several 



cluster expansion [la. 



171]. By making use of the pair approximation 18 . 



19 
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further studies have been concerned with the critical behaviour of the anisotropic spin-1 XXZ 

nn 

Heisenberg ferromagnet with bilinear and biquadratic interactions [18|, [191] , the isotropic spin- 
1 Heisenberg ferromagnet with the bilinear and biquadratic interactions and the single-ion 



anisotropy 



20|, as well as, the anisotropic spin-1 XXZ Heisenberg ferromagnet with an 



antisymmetric Dzyaloshinskii-Moriya interaction [21]. To the best of our knowledge, the 



2 



critical properties of the spin-1 XXZ Heisenberg ferromagnet with the uniaxial single-ion 
anisotropy have not been dealt with in the literature yet. Therefore, the primary goal of 
present work is to examine this model system which represents another eligible candidate 
for displaying an interesting criticality affected by quantum fluctuations. 

The rest of this paper is organized as follows. In Section [ITl we briefly describe the model 
system and basic steps of the variational procedure which gives results equivalent to the 



Oguchi's pair approximation 



22| . Section lllll deals with the most interesting numerical 
results obtained for the ground-state and finite-temperature phase diagrams. Magnetization 
dependences on the temperature, for several values of exchange and single-ion anisotropics, 
are also displayed in Section lllli Finally, some concluding remarks are drawn in Section [IVl 



II. MODEL AND METHOD 

Let us consider the Hamiltonian of the spin-1 quantum Heisenberg model: 

Nq/2 N N 

H = - J J] [A(^f + 5f SJ) + SIS]] - D Y,{Slf -HY, S!, (1) 

where Sf (a = x, y, z) denotes spatial components of the spin-1 operator at the lattice site 
z, the first summation runs over nearest-neighbour pairs on a lattice with a coordination 
number q and the other two summations are carried out over all lattice sites. The first 
term in Hamiltonian ([T]) labels the ferromagnetic XXZ Heisenberg exchange interaction with 
the coupling constant J > 0, A is the exchange anisotropy in this interaction, the parameter 
D stands for the uniaxial single-ion anisotropy and the last term incorporates the effect of 
external magnetic field H . 

The model system described by means of the Hamiltonian ([1]) will be treated within the 
pair approximation formulated as a variational procedure based on the Gibbs-Bogoliubov 



inequahty [231, l2J, |25|, l26| : 



G < Go + - Ho)o. (2) 

Above, G is the Gibbs free energy of the system described by the Hamiltonian ([T]), Gq 
is the Gibbs free energy of a simplified model system given by a trial Hamiltonian Ti^, 
and (. . .)o indicates a canonical ensemble averaging performed within this simplified model 
system. Notice that the choice of the trial Hamiltonian Tig is arbitrary, however, its form 



directly determines an accuracy of the obtained results. If only single-site interaction terms 
are included in the trial Hamiltonian, i.e. single-spin cluster terms are used as the trial 
Hamiltonian, then, one obtains results equivalent to the mean-field approximation. Similarly, 
if a two-spin cluster Hamiltonian is chosen as the trial Hamiltonian, the obtained results 
will be equivalent to the Oguchi's pair approximation, which is superior to the mean-field 
approach. 

In the present work, we shall employ the two-spin cluster approach for the considered 



model system in order to obtain results equivalent to the Oguchi's pair approximation |22 |. 
The two-spin cluster trial Hamiltonian can be written in this compact form: 

N/2 

'Hq = ^^fc, (3) 

k=l 

'^k = ~^[K^kl^k2 + ^kl^k2) + ^kl^k2\ 

-v[isi,r+isi,r]-i{si,+si,), (4) 

where the first summation is carried out over N/2 spin pairs and A, S, rj, and 7 denote 
variational parameters which have obvious physical meaning. It is noteworthy that an 
explicit expression of the variational parameters can be obtained by minimizing the right- 
hand-side of Eq. ([2]), i.e. by obtaining the best estimate of the true Gibbs free energy. 
Following the standard procedure one easily derives: 

A = J, 5 = A, r] = D, -f = {q-l)Jmo + H, (5) 

where tuq = (S'f)o denotes the magnetization per one site of the set of independent spin-1 
dimers described by means of the Hamiltonian T^o- By substituting optimized values of the 
variational parameters ([5]) into the inequality ([2]) one consequently yields the best upper 
estimate of the true Gibbs free energy within the pair-approximation method: 

G = ^G, + ^iq-l)ml. (6) 

Above, Gk labels the Gibbs free energy of the spin-1 Heisenberg dimer given by the Hamilto- 
nian (jll). With the help of Eq. ([6]), one can straightforwardly verify that the magnetization 
of the original model directly equals to the magnetization of the corresponding dimer model, 
i.e. m = (S?) = {S^)q = mo. Of course, similar relations can be established for another 
quantities, as well. 
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To complete solution of the model under investigation, it is further necessary to calculate 
the Gibbs free energy, magnetization and other relevant quantities of the corresponding spin- 
1 dimer model given by the Hamiltonian (jl]). Fortunately, an explicit form of all relevant 
quantities (Gibbs free energy, magnetization, correlation functions, quadrupolar moment) 
can be found for this model system elsewhere [2^]. Referring to these results, the solution of 
the considered model system is formally completed. For the sake of brevity, we just merely 
quote final expressions for the Gibbs free energy Gk and the magnetization mo, both entering 
into Eq. ©: 

Gk = -r'lnZfc, (7) 

Zk = 2exp[/3(A + 2ri)] cosh(2/?7) + 4exp(/?r/) cosh(/37) cosh(/3A(5) 

+ exp[/3(2r/ - A)] +2exp[(3{r] - A/2)] cosh{(3W), (8) 

mo = ^{2 exp[/?(A + 2r/)] sinh(2/37) + 2 exp{pr]) sinh(/?7) cosh(/5A5)}, (9) 
Zd 

where W = ^ {r] - X/2y + 2{X6y, /3 = l/{k^T), k-Q is the Boltzmann's constant, T labels 
the absolute temperature, and the variational parameters A, 5, r/, and 7 take their optimized 
values ([5]). It is quite evident that the magnetization mo must obey the self-consistent 
transcendental Eq. (I9l) (recall that it enters into the variational parameter 7 given by Eq. 
([5])), which might possibly have more than one solution. Accordingly, the stable solution for 
the magnetization mo is the one that minimizes the overall Gibbs free energy (j6]). 

In an absence of the external magnetic field {H = 0), the magnetization tends gradually 
to zero in the vicinity of a continuous (second-order) phase transition from the ordered 
phase (m = mo 7^ 0) towards the disordered phase (m = mo = 0). According to this, the 
magnetization close to the second-order phase transition can be expanded into the series: 

m = am + bm^ + cm^ -(-.... (10) 

Notice that the coefficients a, b, and c depend on the temperature and all parameters 
involved in the model Hamiltonian ([1]). Then, the power expansion of the magnetization m 
can be straightforwardly used to locate second-order transition lines and tricritical points by 
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3^. The 



following the standard procedure described in several previous works [28, 
critical temperatures corresponding to the second-order transitions must obey the condition 
a = 1, 6 < 0, while the tricritical points can be located from the constraint a = 1, 6 = 0, 
and c < 0. Finally, the critical temperatures of discontinuous (first-order) transitions must 



be obtained from a comparison of Gibbs free energy of the lowest energy ordered phase with 
the Gibbs free energy of the disordered phase. 



III. RESULTS AND DISCUSSION 

Before proceeding to a discussion of the most interesting numerical results, let us firstly 
mention that some particular results for the considered model system have already been 



reported on by the present authors elsewhere 33|] . Note that in the former preliminary report 



we have used an alternate approach based on the original Oguchi's pair approximation to 
study a particular case with the coordination number g = 4 corresponding to the square 
and diamond lattices. In the present article, we shall further focus our attention to other 
particular case with the coordination number q = 6, which corresponds to the case of the 
triangular and simple-cubic lattices. A brief comparison with the results obtained previously 
will be made in conclusion. 

Now, let us take a closer look at the ground-state behaviour. A detailed analysis of our 
numerical results shows that the ground-state phase boundary between the ferromagnetically 
ordered and the disordered phases can be allocated with the aid of following condition: 

^ = + (11) 

J 2 g+ 1 ^ ^ 

It is quite obvious from the Eq. ffTTl) that the ground-state phase boundary between the 
ordered and disordered phases shifts to the more positive (weaker) single-ion anisotropies 
when the parameter A is raised from zero. As a matter of fact, the order-disorder transition 
moves towards the weaker single-ion anisotropies for any A 7^ in comparison with the 
result D/J= — g/2 attained in the semi-classical Ising limit (A = 0). This result is taken to 
mean that a destabilization of the ferromagnetic order originates from raising quantum fluc- 
tuations, which work in conjunction with the single- ion anisotropy in the view of destroying 
of the ferromagnetic long-range order at zero temperature. It is worthwhile to remark that 
an appearance of the planar (XY) long-range ordering cannot be definitely ruled out in the 
parameter space with predominant easy-plane interactions {D < and/or A > 1), where 
we have found the disordered phase only. It should be stressed, however, that the present 
form of two-spin cluster mean-field treatment cannot resolve a presence of the ferromagnetic 



long-range order inherent to XY-type models 



34 



35 



36 



371] unlike the conventional Ising- 
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like ferromagnetic long-range order with only one non-zero component of the spontaneous 
magnetization. 

Next, let us turn our attention to the finite-temperature phase diagram, which is shown in 
Fig. 1 in the reduced units d = D / J and t = k^T/ J for the simple-cubic (triangular) lattice 
and different values of the exchange anisotropy A. In this figure, the solid and dashed lines 
represent second- and first-order phase transitions between the ferromagnetic and paramag- 
netic phases, respectively, while the black circles denote positions of the tricritical points. It 
is quite obvious from this figure that the considered model system exhibits the highest values 
of critical temperature in the Ising limit (A = 0). The gradual increase of the exchange 




FIG. 1: The phase diagram of the spin-1 Heisenberg model for the simple-cubic (triangular) lattice 
and several values of the exchange anisotropy A. The solid and dashed lines represent second- 
and first-order phase transitions, respectively. The black circles denote positions of the tricritical 
points. 



anisotropy A reduces the transition temperature as a result of raising quantum fluctuations. 
It is worthwhile to remark that all the lines of second-order phase transitions, for arbitrary 
but finite A, have the same asymptotic behaviour in the limit d —>■ oo. Actually, the critical 
temperature of the continuous transitions does not depend on the exchange anisotropy in this 
limiting case and it is equal to t* = 5.847. Moreover, it should be also mentioned that our 
approach yields for the Ising case without the single-ion anisotropy (A = 0, = 0) the crit- 
ical temperature tc = 3.922, which is consistent with the result of other pair-approximation 



methods 



2l| and is simultaneously superior to the result tc = 4.0 obtained from the stan- 
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FIG. 2: The phase diagram of the spin-1 Heisenberg model for g = 6 and A = 0.0, 1.0, 2.0, 3.0, and 
4.0. The solid and dashed lines represent second- and first-order phase transitions, respectively. 
The black circles denote positions of the tricritical points. The dot-and-dash line represents the 
location of tricritical points in dependence on the exchange anisotropy A. 



dard mean- field approximation 38|]. In addition, it can be clearly seen from Fig. 1 that 
the transition temperature of the continuous phase transition monotonically decreases by 
decreasing the single-ion anisotropy d until the tricritical point (TCP) is reached. Further 
decrease of the anisotropy parameter d changes the second-order phase transitions towards 
the first-order ones. It should be realized, nevertheless, that the first-order phase transitions 
occur merely in a narrow region of single-ion anisotropics close to the boundary value (fTTl) 
at which both completely ordered phases with m = ±1 have the identical energy (coexist to- 
gether) with the disordered phase with m = and one asymptotically reaches the first-order 
phase transition between them in the zero-temperature limit. An origin of discontinuous 
phase transitions could be therefore related to the fact that the ordered and disordered 
phases have very close energies near the boundary single-ion anisotropy (ITT!) (the former 
ones are being slightly lower in energy) and the small temperature change might possibly 
induce a phase coexistence (energy equivalence) between them, what consequently leads to 
the discontinuous phase transition. In Fig. 2, we depict more clearly the position of TCPs 
in dependence on the single-ion and exchange anisotropics by the dot-and-dash line in order 
to clarify how the type of phase transition changes with the anisotropy parameters. As one 
can see from this figure, the d-coordinate of TCPs {dt) shifts to more positive values upon 
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FIG. 3: The temperature dependence of the magnetization m for the isotropic spin-1 Heisenberg 
model (A = 1.0) on the simple-cubic (triangular) lattice, when the value of the single-ion anisotropy 
parameter d changes. The dashed lines represent the discontinuities of the magnetization at the 
first-order phase transitions. 

strengthening of A, while the ^-coordinate of TCPs behaves as a non-monotonic function of 
the exchange anisotropy A with a minimum at A^j„ = 3.459. 

To illustrate the effect of uniaxial single-ion anisotropy on the phase transitions, the 
thermal variation of the magnetization m is shown in Fig. 3 for the case of isotropic spin-1 
Heisenberg model (A — 1.0) and several values of d. It can be clearly seen from this figure 
that the reduction of the single-ion anisotropy causes lowering of the critical temperature 
tc- Furthermore, it is also evident that the magnetization varies smoothly to zero for d — 
0.0, —2.0, and —2.6 until the temperature reaches its critical value. This behaviour of 
magnetization, which is typical for the second-order (continuous) phase transitions, persists 
until d > dt {dt = —2.656 for A = 1.0 and q = 6). On the other hand, the magnetization 
jumps discontinuously to zero for d < dt (e.g. see the curves for d = —2.7 and —2.8), what 
is characteristic feature of the first-order (discontinuous) phase transitions. As one can see, 
this discontinuity in the magnetization increases rather abruptly as the single-ion anisotropy 
moves to more negative values with respect to the dt value. Finally, it should be pointed 
out that the similar variations of magnetization curves occur for any value of the exchange 
anisotropy A. 
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IV. CONCLUSION 



In the present paper, the phase diagram of the anisotropic spin-1 XXZ Heisenberg model 
with the uniaxial single-ion anisotropy is examined within the variational procedure based 
on the Gibbs-Bogoliubov inequality, which gives results equivalent to the Oguchi's pair ap- 
proximation [22]. A comparison between the results obtained in the present study and those 
attained within the standard Oguchi approximation actually implies an equivalence between 
both the methods. The most important benefit of using the variational approach based on 
the Gibbs-Bogoliubov inequality is that in this way adapted method enables obtaining of 
all thermodynamic quantities in a self-consistent manner and moreover, it is also well suited 
to discern the continuous phase transitions from the discontinuous ones by distinguishing of 
the stable, metastable and unstable solutions inherent to the approximation used. 

In the spirit of the applied pair-approximation method we have demonstrated that the 
single-ion anisotropy as well as the exchange anisotropy have a significant influence on the 
critical behaviour and both these anisotropy parameters can cause a tricritical phenomenon, 
i.e. the change of the continuous phase transition to the discontinuous one. Our results 
can serve in evidence that the tricritical phenomenon may occur in the investigated model 
system if at least one of the anisotropy parameters provides a sufficiently strong source of 
the easy-plane anisotropy. Note furthermore that the obtained results are rather general in 
that they are qualitatively independent of the lattice coordination number. The comparison 
between the results to be presented in this work with those reported on previously for other 



particular case [33| actually implies that the model under investigation shows qualitatively 



the same features irrespective of the lattice coordination number. 
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